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Summary. In these lectures I introduce aspects of current theory used to inter- 
pret the preUminary data on ultra-relativistic nuclear collisions at RHIC energies 
in terms of the physical properties of QCD matter at extreme densities. Topics 
covered include: What are the physics questions at SPS and RHIC?, Geometrical 
vs Dynamical features of A+A, the interplay of computable Hard pQCD vs phe- 
nomenological Soft dynamics, Baryon number transport and Junctions, How can 
we compute and get experimental control over the initial conditions?, how to rec- 
oncile apparent hydrodynamic behavior with partonic/hadronic transport theory. I 
use the preliminary RHIC data available up to June 1,2001 to illustrate these top- 
ics. Most technical details are deferred to the literature. However, since the main 
new observable at RHIC relative to SPS is jet quenching, I elaborate more on this 
"tomographic" probe of ultra-dense matter. The possible discovery of jet quenching 
at RHIC by STAR and PHENIX is highhghted. 



1 Introduction 

Finally, after 20 years of preparation |Q, a new chapter in nuclear/particle 
physics commenced on June 12, 2000 with the measurement of the first 
Au + Au collisions at -/s — 56 AGeV (GeV per nucleon pair) in the Rel- 
ativistic Heavy Ion Collider (RHIC) at the Brookhaven National Lab (BNL). 
Soon thereafter collisions at ^/s = 130 AGeV were also measured. The first 
results were reported at Quark Matter 2001 Q from the four major experi- 
ments, STAR||, PHENIxg, PHOBOSg, and BRAHMS[||. A small army 
of ^ 1000 experimentalists measured the flavor, rapidity and transverse mo- 
mentum distributions of the approximately 4000 charged particles produced 
in each central (head on) collision at 130 AGeV. In the summer of 2001, it is 
anticipated that RHIC will reach its design energy, and p + p and Au + Au 
collisions at ^/s — 200 AGeV, will come under experimental scrutiny. 

These lectures provide a very condensed introduction to current theoreti- 
cal work aimed to provide a consistent interpretation of observables measured 
in such reactions in terms of the properties of dense QCD matter. The color 
slides of the original lectures can be found on my WWW site 0. These lec- 
tures note are designed to supplement those slides and update them with the 
preliminary RHIC data available as of June 1, 2001. 

The theoretical work on the new physics that may exist in QCD matter 
at extreme densities began in the mid 1970's with the realization that the 
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asymptotic freedom property of QCD implies the existence of a new phase 
of strongly interacting matter called the Quark-Gluon Plasma (QGP) 
[ p^ . Unlike familiar nuclear or more generally hadronic matter consisting of 
composite "elementary" particles {tt, K, p,p, A, A, ■ ■ ■) in which quarks and 
gluons are permanently confined, the QGP phase at very high temperature 
and/or baryon chemical potential {T,fiB ^ Aqcd ~ 200 MeV = 1/fm) is 
one where the interactions between quarks and gluons become relatively weak 
and short range 

V(r) ^ , Qfc = — oc - — — . 

^ ' An r ' 47r log(ror/^B) 

The color electric (Debye) screening mass fioiT, Hb) increases linearly with T 
OT hb modulated by a slowly varying factor of the running coupling, g{T, ^b) 
(see lectures of Rebhan). The thermodynamic properties of this dcconfincd 
and chirally symmetric (~massless) phase of matter are thus expected in 
perturbation theory to reduce approximately to an ideal Stefan-Boltzmann 
gas of quarks and gluons. For the Standard Model with 3 colors and Nf 
flavors of "light" quarks relative to T, (SUdS) (g) SUf{Nf)), the Stefan 
Boltzmann constant for /is = is 

i^ss = ^ = ^ = (2. X 8e + ^ X 2, X 2,+,- X 3, X iV/)^(l + 0(5')) - 12 

taking the helicity, color, flavor, and antiquark degrees of freedom into ac- 
count. In reality the severe infrared singularities of perturbative QCD (pQCD) 
lead to large non-perturbative corrections to the ideal gas equation of state for 
all temperatures and chemical potentials accessible experimentally even be- 
yond the future Large Hadronic Collider. Only numerical lattice QCD (IQCD) 
methods (see lectures of F. Karsch) can provide reliable predictions for 
the thermodynamic properties of the QGP phase of matter. Effective models 
and resumed many-body techniques (see Alfors, T. Rebhan, J. P. Blaizot, and 
E. Shuryak) are however needed to interpret the IQCD "data" and provide 
physical insight especially at finite chemical potential. However, it is sobering 
to recall that pQCD thermodynamic expansion of the pressure in powers of 
g shows no sign of convergence even before the Linde infrared catastro- 
phe at 0{g^), and non-perturbative corrections to the pQCD Debye mass, 
Hd, remain about a factor of 3 up to T ^ 200 GeV[|l^. The full the- 
oretical understanding of the structure of the non-abelian plasma phase of 
QCD therefore remains a fundamental open problem in physics because it in- 
volves strongly correlated, nonperturbative and possibly turbulent dynamical 
features p5|. 

One of the essential and intriguing aspects of the QCD many-body prob- 
lem is that the physical vacuum is an extraordinarily complex coherent many 
body medium. The gluon and quark condensates lower the energy density of 
the physical vacuum by and amount B ^ ^qcd ~ "^^^ MeV/fm'^. Drilling a 



Lectures on the Theory of High Energy A+A at RHIC 



3 



perturbative vacuum bubble of volume V in this condensate costs an enor- 
mous energy BV. The QGP, if formed in V, must counteract the physical 
vacuum pressure, B to prevent it from collapsing. This is only possible when 
the temperature exceeds « {3B/KsbV^'^ 150 MeV. 

The driving force behind the experimental effort at CERN and BNL over 
the past 20 years has been to try to create the extreme conditions necessary 
to produce and diagnose this new form of matter in the laboratory. Over 
the past 15 years experiments at the AGS/BNL {^/s = 5 AGeV) and the 
SPS/CERN {y/s = 20 AGeV) have searched systematically via a very large 
set of observables for evidence of the QGP phase (see lectures of C. Lourenco). 
In these lectures I will focus on the most recent developments in that search 
that has just begun at RHIC. I must emphasize from the onset that most 
of the data shown here are of PRELIMINARY nature and could change as 
better control over the systematic errors is achieved in the next few year. 
Nevertheless, the new data are so exciting that it is worth trying a first pass 
to give an overview and possible interpretation. 

2 Geometry and Dynamics in A+A 

The main obstacle in interpreting data on collisions of finite nuclei (at any en- 
ergy) is that the matter created undergoes quantum (perhaps semi-classical) 
many-body dynamics that may be approximated by thermodynamics only 
over a limited (low px) kinematic range. Experimentalists do not have the 
luxury of lattice or perturbative QCD theorists of tapping into the infinite 
gedanken volume or reservoir with a fixed temperature and pressure. Nuclear 
collisions produce dense matter in a highly dynamical environment, and the 
matter produced expands anisotropically near the speed of light. It is far 
from clear whether local thermal and chemical equilibrium concepts apply, 
and even so, over what domain of the 8 dimensional {x^ ,p^^) phase space can 
they be used. 

Before a collision, the partons of the two colliding nuclei are locked into 
a coherent field configuration. The dense virtual cloud of gluons and quarks 
may be described in the colinearly factorized QCD approximation by A times 
the known structure functions, fa/pix,Q), of nucleons when the resolution 
scale is high enough Q > 1 — 2 GeV. However, many body initial state inter- 
actions could lead to strong modifications of this naive parton picture (see 
lectures of L. McLerran). The nuclear QCD fields continue to interact after 
the nuclear valence quark pancakes pass through each other. The interaction 
spans a space-time hyperbola over a proper time ^Jt"^ — ~ 30 fm/c = 10~^^ 
sec. Then a "miracle" happens! The field quanta hadronize in way that is 
unfortunately not well understood. The dense final hadronic debris can fur- 
ther interact as it expands toward the detector elements. From CERES/SPS 
data|p^ there is evidence that the in-medium mass-width (spectral function) 
of vector mesons may change drastically Q. 
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Hard probes (jets, leptons, photons, heavy quarks) are of special inter- 
est because they provide "tomographic" tools with which can map out this 
evolution experimentally. Hard probes are effective "external" tomographic 
probes because they are produced with a pQCD computable initial distri- 
bution on a much shorter time scale, l/mj_, then the plasma formation 
time, ~ 1/3T. Modification of their known initial distributions therefore pro- 
vides information on the medium through which they propagate in analogy 
to conventional X-ray or positron tomography used in medicine |p^ . The pri- 
mary advantage of RHIC over lower energy machines (AGS, SPS) is that 
hard pQCD probes are produced at RHIC orders of magnitude more abun- 
dantly over a significantly larger kinematic range. This greatly improves their 
tomographic resolution power. 

Figure 1 shows the rapid growth of high pr Au + Au +X predicted 
by pQCD from SPS to RHIC and LHC. As discussed for example in it 

Au+Au (b<3) -> 7t° s = 20, 200 , 5500 AGeV 
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Fig. 1. Invariant distribution of hard pQCD produced ^° in central Au -\- Au colli- 
sions as a function of cm. energy via eq.(Q,^ with CTEQ5M structure functions, 
K=2 factor, scale Q = Pt/'2, and multiplied by nuclear overlap Tab ~ 24/mb. The 
dashed curve shows the contribution from gluon jet fragmentation only. 




is useful to decompose the nuclear geometry dependence of invariant hadron 
distributions produced inA + B—th + Xat impact parameter b into a phe- 
nomenological "soft" and perturbative QCD calculable "hard" components 
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as 



dNAB jh) _ f,. dN,oftih) 1 dahardih) 

d-^p aya^PT <„ dyd-^ pt 

where NpartO^) is the number of nucleon participants and Ncou(h) = a^^TAsO^) 
is the number of binary NN coUisions at impact parameter b. The nuclear 
geometry of hard collisions is expressed in terms of the Glauber profile density 
per unity area TA_B(b) = / d^r T^(r)TB(r — b) where TA(r) = J dz Pa{tc, z) 
(see Fig. 2). The hard part scales with the number of binary collisions cx A^^^ 
because their probability is small built up from all possible independent par- 
ton scattering processes. The soft part scales with only Npart oc because 
their probability is large and therefore "shadowed" . 



A+A Collision Geometry Experimental Handles 



Participant Nucleon Number 




Fig. 2. Illustration of key aspects of the relation between the geometry of nuclear 
collisions and the participant and collision number at a fixed impact parameter. 
The observables (see[^Q|5|Q) used to constrain the geometry experimentally are 
also illustrated. 



The (textbook) computable lowest order pQCD differential cross section 
for inclusive p + p —t h + X invariant cross section is given by 



^^^^pL =KY, JdXadXbfa/pi00a,Ql)h/p{xb,Ql)^iab ^ Cd)- 
abed 



(2) 
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where Xa ~ Pa/PA,Xb = Pb/Ps arc the initial momentum fractions carried 
by the interacting partons, Zc — Ph/Pc is the momentum fraction carried by 
the final observable hadron, fa/p{xa,Qa) is the proton structure function 
for parton of fiavor a, and Dfi/c(zc,Qc) is the fragmentation function for 
the parton of flavor c into h. The UAl data on pp hadron production with 
Pt > 2 GeV can be well reproduced with this pQCD model expression. 

The soft {pt < pq ^ 2 GeV/c) nonperturbative contribution to the 
hadron yields can only be modeled phenomenological. The Dual Parton Model 
[^2|2^ and the LUND string model ||2j,^ are the most extensive and success- 
ful low Pt multiparticle phenomenologies. The basic pQCD matrix elements 
have been encoded into a Monte Carlo code, PYTHIA |26|. A variant of soft 
string phenomenology tuned to pp,pp data, with the hard part taken from 
PYTHIA, a hadronization scheme taken from the LUND JETSET hadroniza- 
tion, and a eikonal nuclear multiple collision geometry were combined into the 
Monte Carlo A+B collision generator in HIJING (2^. HIJING has been used 
over the past decade to predict many observables at RHIC po| , ^7|j2^1 . The 
separate soft and hard components in HIJING with a fixed A, indepen- 
dent scale po = 2 GeV/c are illustrated in Fig. 3. Hard gluons in the LUND 

TWO COMPONENT MODELS FOR A+A 

1) Np3rt(b) Soft Beam Jets (LUND Strings) 

— tud 
ud * • u 



2) N,,„(b) Hard pQCD 

Pj>Po=2 GeV/c (PYTHIA) 



3) = ^+2 Hard pQCD Hadronized in Strings 

(HIJING) u«--<Tr"~^ud 
ud ^^bZ^^ u 



Fig. 3. Illustration of Hard-f Soft two component models of A-l-A reactions like 
HIJING @. 

hadronization scheme are represented by kinks in the strings between valence 
quarks and diquarks of the Npart (b) interacting baryons in A -\- B collisions 
at impact parameter b. In this class of models no final state interactions are 
taken into account. 

The physics in A + A reactions that must be understood in order to be 
able to interpret observables in terms of the properties of dense QCD matter 
requires extending the above class of event generator models to include 
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1. Initial Conditions: The formation physics responsible for creating an in- 
coherent gas of giuons and quarks from the initial virtual nuclear fields, 

2. Parton Transport: The {x'^,p'^) phase space evolution of that parton gas 
toward equilibrium, 

3. Hadronization: The dynamical mechanisms that convert the parton de- 
grees of freedom in the confining physical vacuum into the observable 
hadronic ones, 

4. Hadron Transport: The final state interactions of the expanding dense 
hadronic matter prior to "freeze-out" . 

Each problem is fascinating in its own right but only bits and pieces are un- 
derstood or phenomenologically mapped out up to now. There exists unfortu- 
nately no complete computable dynamical theory (like Magneto-Hydrodynamics 
for QED plasmas) that consistently takes into account all four elements. QCD 
is believed to be THE theory, but it is still not computable except at high px 
where perturbative or classical methods may apply. There exists instead sev- 
eral diff'erent dynamical "scenarios" to describe A + A that attempt to patch 
together different approximation techniques and phenomenological models to 
address all the physics issues in turn. 

Two generic approaches to A + A can be classified by whether the Initial 
Conditions are (1) computed (via pQCD or classical Yang Mills (cYM)) and 
subsequent evolution followed by a dynamical scenario for 2-4, or (2) the 
initial conditions are fit by extrapolating final observables backwards via a 
suitable dynamical scenario. At lower energies (AGS,SPS) only the second 
approach is available since the momentum scales are simply too low to apply 
either pQCD or cYM. At collider energies RHIC and beyond, the copious 
production of mini -jets |2^,|^,||j35|j3|l with pT > po ^ 2GeV shown in Fig.l 



makes it possible for the first time to pursue the first approach via pQCD 
Eq.(|2|). At very high energies classical Yang Mills theory ||^,|8|j39|j4^,|l|,|2| 



provides a general method to compute the Formation Physics which reduces 
to pQCD at high pT- Whether RHIC or LHC energies are high enough is an 
open question. 

The second approach, trying to "fit" the initial conditions by extrapo- 
lating the final distributions backwards with a suitable dynamical model has 
been traditionally based on relativistic hydrodynamics p^ , p^ , [45[|4^ . The ap- 
proximate longitudinal boost invariant boundary conditions at ultra-relativistic 
energies simplify hydrodynamic equations greatly as pointed out by Bjorken|^^ . 
For /is = the hydrodynamic equations are, 

d^T^'" = , T^'^ix) = u^'u''{e + P) - g^'^P , 

where e(x), P{x) are the proper energy density and pressure and u^{x) is the 
four velocity field of the fluid. The central assumption is that thermal and 
chemical equilibrium are maintained locally in spite of the possible large gra- 
dients in the fluid variables. The great advantage of hydrodynamics is that 
it provides a covariant dynamics depending only on the equation of state 
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P{T{x)) that is directly related to the IQCD predictions. When a specific 
space-time freeze-out hypersurface is assumed together with the assumption, 
the Cooper-Frye prescription p9| , p0[j5l| , |5^ , p^ , the computed four fluid veloc- 
ity field can be used to predict the final anisotropic fiow pattern of hadrons. 
Since this process is assumed to be reversible, the final distributions together 
with an assumed freeze-out hypersurface can be used to compute the ini- 
tial conditions on any desired initial hypersurface. The disadvantage of this 
approach as emphasized in ||55|] is that both the initial and final freeze-out 
hypersurfaces must be guessed. Also finite mean free path physics is outside 
the scope of ideal hydrodynamics, and transport theory solutions|5^ do not 
support "sharp" freeze-out hypersurfaces. Thus the inversion of data in this 
way to deduce the initial conditions is not unique. The neglect of dissipa- 
tive effects such a viscocity also makes it impossible to relate central A + A 
to peripheral and light ion data, especially p + p. Finally, the assumption 
of homogeneous or slowly spatially varying initial conditions is questionable 
because of copious mini-jet production In spite of all the above theoreti- 
cal problems, initial conditions for RHIC have been successfully constructed 
that lead via ideal hydrodynamics and idealized Cooper-Frye freeze-out to 
distributions that reproduce amazingly well many of the low pr observables 



at RHIC|53,|58|,|9|J60| (see next section). 

In order to bring the freeze-out assumption under better theoretical con- 
trol covariant, nonequilibrium transport theory]^ must be solved. Until re- 
cently, only simplified 1-l-lD Bjorken transport theory was soluble in the 
linearized relaxation time approximation (see ]62| and refs therein). This is 
due to the great numerical complexity of the 3+lD nonlinear Boltzmann 
equations p5| , |63t : 

P'tdji - jjjifzfA- flf2)Wi2^^i5\pi+P2-pz~PA)+ S{x,pi). (3) 
234 



where W is the square of the 2^2 scattering matrix element, the integrals 
are shorthands for / = / ^2^ye- ' "^here g is the number of internal degrees 

i 

of freedom, while fj = /(a;,pj) is the parton phase space distribution. The 
initial conditions are specified by a source function 5(a;,p) that corresponds 
to the assumed initial conditions. 



Yang Pang's parton subdivision technique [p 3 64 and the speed of cur- 
rent workstations have finally made it possible to solve eq.(^ numerically. 
(Codes can be obtained from the OSCAR Web site|65j). The solutions | |66| , p7t 
prove that elastic parton scattering with pQCD rates is insufficient at RHIC 
to keep the plasma in local equilibrium due to the extreme rapid longitudinal 
"Hubble" expansion of the system|6^. Unfortunately, there exists no prac- 
tical algorithm at this time to solve the more nonlinear inelastic transport 
equations involving gg ng processes. Therefore, if hydrodynamics applies 
to A + A at RHIC, then most likely strong nonpcrturbative mechanisms must 
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be assumed to exist on faith or hypothesis (see E. Shuryak lectures). This is 
an important open theoretical problem. 

I would also like to call attention to a new class of hydrodynamic model p8| 
that side-step the final freeze-out problem by assuming that local equilibrium 
is maintained only up to an intermediate hyper-surface, just after hadroniza- 
tion on a T = Tc — e isotherm. Using that intermediate freeze-out as the 
initial conditions of a hadronic transport theory, the subsequent evolution of 
the hadronic gas toward a dynamical freeze-out is then determined by known 
hadronic cross sections via URQMD [|69y70[] . 

3 Preliminary Results from RHIC 
3.1 Global Constraints on Initial Conditions 

The first results from RHIC from PHOBOS||j, shown in Fig.|[ demonstrate 
that the energy dependence of the scaled charged particle (pseudo)rapidity 
density, dNch/drj/Npart, is different from p + p and p + p systematics. 
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Fig. 5. Normalized pseudorapidity 
density (PHOBOS and PHENIX) g] 

diVch/d?7||^l<i/(0.5 X Npart) as a 
function of the number of partici- 
pants. Predictions based on HIJING 
(thick solid) and EKRT|||| (thin) are 
shown. 



Fig. 4. Measured pseudorapidity 
density normalized per participant 
pair for central Au+Au collisions 
(PH0B0S§j2|]). Systematic errors are 
shown as shaded areas. Also shown are 
results of Pb-hPb data (CERN SPS), 
HIJING 1^ simulations and a parame- 
terization of pp data. 

Approximately 50% more particles are produced at mid rapidity per par- 
ticipating baryon in central Au+Au collisions then in p+p at the same energy 
per baryon. The curve shows that the two component HIJING model pre- 
dicted well this result. However, as shown in Fig||, another model EKRTJsgI, 
was also found to predict the same multiplicity as HIJING for central col- 
lisions. In ref. ||2^, we proposed that the centrality dependence of this ob- 
servable could differentiate between these competing models of the initial 
conditions. The new data of PHENIX and PHOBOS verified this predic- 
tion. While neither model accounts quantitatively for the data, but the two 
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component HIJING model with its combined and A^^"^ dependence bet- 
ter describes the rate of increase of the scaled multiplicity with participant 
number. The observed increase of the scaled multiplicity with energy relative 
to p + P and with participant number dependence is supports the prediction 
of copious mini-jet production at RHIC. This is one of the necessary, though 
insufficient, conditions to form a dense gluon plasma in A -I- ^. 

The difference between HIJING and EKRT is that in the latter it is 
assume that all the produced entropy (multiplicity) arises at RHIC energies 
from hard pQCD processes. EKRT assume that there is no significant soft 
component, i.e. dNsoft ^ dNhard in eq.(|l]). However, the hard component is 
cutoff at scale pa that is allowed to vary with both A and hence Npart and 
with energy -^s based on the following assumption: independent and hence 
TAB{b) proportional number of gluons with px > Po are produced only in 
"resolvable" domains of finite area tt/pq. There are PqR^ such domains in 
the transverse plane in a central nuclear collision. This so called "final state 
saturation" model is then specified by 

^ - r - pviR^ ■ (4) 



For (3=1 assumed in EKRT, the solution for the saturation scale is po {y/s, A) e 
Psat « 0.2A° i3(^s)0 i9_ This predicts dNg/dy cx A°-^^ in spite of the appar- 
ent proportionality of hard processes to A'^/^ . The flat (dNg/dy) /Npart A~° 
independence of the scaled multiplicity is a general feature of saturating QCD 
models of the initial conditions (see also lectures of Mclerran). Such a flat 
behavior is however ruled out by the present data at RHIC. 

An alternate (so-called initial state saturation) model was proposed by 
KN|Q based on the nonlinear QCD evolution equations of In this model 
of nuclear initial conditions, the number of liberated gluons is proportional to 
the number of virtual gluons participating in the reaction on a scale pq. The 
produced number is then taken to fNpartxG{x,pQ) in terms of the nucleon 
gluon structure function, where / ~ 1.2 is a factor on the order of unity. Since 
the interaction probability is proportional to the running coupling as{po), the 
initial state saturation condition is defined by 

^ = fNpartxG{x,po) = f-^^—-plR'. (5) 
dy 3n''as(po) 

The main difference between initial and final state saturation models is there- 
fore due to the logarithmic dependence on pq introduced by the running cou- 
phng. In a simple ansatz was assumed for xG{x, Q) oc logQ/A based on 
the linear (DGLAP) evolution equations. With this ansatz KN predicted a 
participant dependence surprisingly close to the observed data in Fig.^. 

However, the x independent ansatz of KN used for xG{x, Q) for the scale 
Q ^ 1 GeV/c is a guess that cannot be supported by the available ep HERA 
data. At small x ~ 0.01 and low Q ~ 1 the pQCD factorization analysis of 
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deep inelastic e+p reactions breaks down and xG acquires a 100% systematic 
error bar as shown in Fig.^. Initial state saturation is a theoretically sound 
model only at very high energies or nuclei with A ^ 200, when Q > 2 GeV/c 
and the errors based on pQCD analysis become manageable. 



Initial State Saturation at E^„=130 AGeV 
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Fig. 6. Graphical solution to the initial state saturation eq.( |^ for different impact 
parameters at 130 AGeV using the xG gluon structure functions from ZEUS and 
CTEQ. The energy dependence enters through the dependence on Xsat = "iQsatl \fs. 
The participant number dependence of Qaat follows from the intersection of the 
parabolic curves with xG. Unfortunately Au nuclei are too small, and the solutions 
in the Q ~ 1 GeV region are completely unreliable. 

While it is premature to conclude which approach is least wrong (see 
also [^S)), in my opinion, it appears that the global multiplicity data and 
its centrality dependence can be used as indicators that the initial gluon 
rapidity density at RHIC is between HIJING's 200 and EKRT's 1000. The 
corresponding gluon density, Pg{T) = dNg/dy/iTTiE?) is thus ^ 10 — SO/fm'^ 
at the corresponding formation time 1/po ~ 0.1 — 0.2 fm/c. Thus RHIC may 
have indeed created the densest gluon plasma ever in the laboratory. As I 
emphasize in a later section, fortunately there are many other observables, 
especially jet quenching, that provide independent checks of this possibility. 

It is important to emphasize that similar results for the multiplicity in 
central collisions in HIJING and EKRT are purely coincidental because the 
models differ by a factor of five on the initial gluon density. This is compen- 
sated for by the underlying very different hadronization schemes assumed. 
HIJING creates a large fraction of the observed hadrons at RHIC through its 
soft string fragmentation scheme, while EKRT assume that entropy conser- 
vation implies that cLNt^ « dNg. The lack of a detailed hadronization theory 
can only be overcome phenomenologically by testing experimentally all the 
ramifications any particular model. 
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Another observable that was suggested in ||22l to help differentiate mod- 
els of initial conditions is the shape and scaling of the whole rapidity dis- 
tribution (see Fig.R). It is seen that HIJING predicts a somewhat narrower 




Fig. 7. (a). Measured dNch/dr]/{{Npart)/2) for {A'j,art) = 102 (circles), 216 
(squares) and 354 (diamonds) by PHOBOS|3o|. (b). Same as (a) from HI- 
JING. 

and stronger centrality dependence than observed by PHOBOS. This may be 
related to the baryon stopping power at RHIC. Unfortunately no predictions 
are available for either the initial or final saturation models on the predicted 
shape of the rapidity distribution. This observable is especially sensitive in 
those models to the x dependence of the saturation criteria. 

The total integrate charge particle multiplicity is shown in Fig. ^. RHIC 
has produced about 4000 charged particles in Au+Au at 130 AGeV. The non- 
linear enhancement near central collisions is interpreted in terms of HIJING 
as due to the onset of the mini-jet component. 



3.2 Global Barometric Observable Et/Nch 

An important global barometric measure of the internal pressure in the ultra- 
dense matter produced is the average transverse energy per charged particle. 
PHENIX data are shown in Fig.|| compared to WA98 data from CERN. 
What is most amazing is that Er/Nch ~ O.SGeT^ almost independent of ^/s 
from 20 to 130 AGeV and independent of centrality! HIJING predicts that 
it should rise from 0.8 to 0.9 GeV from CERN to SPS due to the enhanced 
mini-jet activity at RHIC. The EKRT initial state saturation model predicts 
a growth of this quantity in the initial state by about a factor of 3. The 
reason that EKRT remains viable after these data is that the assumed entropy 
conservation implies that a large amount of pdV work due to longitudinal 
expansion is performed by the plasma. In 1-l-lD hydrodynamics the energy 
per particle e/p w 2.7T decrease as the system expands and cools T ^ 
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400 



Fig. 8. PHOBOS total charged particle multiplicity vs nucleon participant 
numberlpsl. 




100 150 200 250 300 350 400 
Npart 

Fig. 9. Preliminary PHENIX data and WA98 data on dEr/dNch as a function of 
participant number[^. This barometric observable appears to be independent of 
= 20 - 200 and centrality! 



if the freeze-out is assumed to occur at all energies and impact parameters 
in ^ + j4 on a fixed decoupling isotherm, then the energy per particle will 
always be the same. At RHIC this global transverse energy loss from the 
initial state is predicted to be about a factor 3. The theoretical problem of 
justifying hydrodynamics and the freeze-out prescription itself discussed in 
the previous section comes back to haunt us here|^. The observed NULL 
effect in Ex/Nch is very interesting because it is so difficult to obtain in any 
transport theory with finite pQCD relaxation rates. 
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3.3 Discovery of Jet Quenching 

One of the predicted |2|j7^,Q signatures of dense matter formation is the 
suppression of jets and their high pT hadronic debris due to energy loss of 
the jet in the medium. However, the search for this effect at SPS by WA98 
yielded the opposite result as shown in FigJlG. Even a modest dEdx = 0.2 
GeV/fm is completely ruled out by the data[|81[. The problem is that at lower 
energies, multiple initial state elastic scattering leads to a random walk in 
transverse momentum. This enhances the pt of the scattered partons so that 
(Pt) = + A^/^Sp^. This so-called Cronin effect has been well studied in 
p + A reactions up to 800 GeV. At lower energies the very steep fall of the 
high Pt tail makes the distribution extremely sensitive to this modest px 
enhancement. When convoluted through two nuclei, Wang predicted[^ that 
the Cronin enhancement at SPS in Pb+Pb should be a factor of two as verified 
in Fig.|ri|. What is plotted there is the ratio of the observed invariant cross 
section to the scaled binary collision number, Ncou{b), scaled invariant cross 
section in p+p. Unity corresponds to naive superposition of NcoU independent 
elementary p + p hard processes in the absence of any nuclear effects. The 
ratio starts below 1 since the low pt distribution grows only with the number 
of participants (divided by two) and Npart{0) f^NcoiiiO) « 0.15. 

In stark contrast to the SPS enhancement of high px pions, a factor of 
two or more suppression of > 2 GeV hadrons was reported by STAR|^,U| 
and PHENIX|8|. 




m^-rriQ (GeV) 

Fig. 10. Single-inclusive vr" spectra in central S+S at Euh = 200 GeV and Pb+Pb 
collisions at Siab = 158 GeV. The solid lines are pQCD calculations (Wang[^) with 
initial-^T broadening and dashed lines are without. The S + S data are from WA80 
and Pb + Pb data are from WA98 The dot-dashed line is obtained from the solid 
line for Pb + Pb by shifting pr by 0.2 GeV/c. 

Fig.^ shows that for px < 2 GeV a similar trend of increase due to the 
gradual transfer from participant to binary scaling is taking place as at SPS, 



Lectures on the Theory of High Energy A+A at RHIC 



15 



Q. 
+ 

A 



^1 
\ 

m 
+ 

< 



E„=158AGeVPb+Pb 

A n" 10% central 
t h' 5% central (NA49) 
i n 8% central Pb+Au (CERES) 




P,(GeV/c) 

Fig. 11. The nuclear modification 
factor for hadron spectra in central 
Pb + Pb collisions at the CERN- 
SPS exceeds unity at high pr due to 
the Cronin effect. The solid line is a 
pQCD calculation by WangQ. 
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Fig. 12. The nuclear modification factor 
for charged hadrons in central Au + Au 
RHIC from STAR||4). In contrast to 
SPS, the high pr charged hadrons are 
suppressed. 



but for pT > 2 GeV the ratio for charged particles tt^ + + starts to 
drop again and reaches ~ 0.5 at 4 GeV/c. 



The PHENIX data 
for identified ttO in Fig 



shows an even more dramatic quenching pattern 
.13. In this experiment, it was further verified that 
"peripheral" colHsions are not quenched while central ones are. Fig.|l^ shows 
that the suppression factor may reach a factor of three at 3 GeV/c. In this 
plot the ratio is not relative to pp data extrapolated to 130 GeV, but to 
"peripheral" collisions where the average number of participants and binary 
collisions is only « 20. In contrast Npart ~ 360 and Ncoii ~ 857 for the central 
collisions. It must be emphasized that current systematic errors are still much 
larger than statistical, but it is clear that the combined information from two 
independent experiments in Figs.( p"2|jl4|) imply that something new has been 
discovered in A + A collisions at RHIC. I believe that this is the predicted 
jet quenching as discussed in the next lecture. 

The reason that this discovery is perhaps even more exciting than the 
famous J/ip suppression effect discovered by NA50 p6| , ^ at the SPS is that 
J/ijj suppression was also seen in p + A. The cold nuclear suppression mech- 
anism in p + A is called "normal" . The enhance suppression in Pb + Pb is 
"anomalous" because it is more than if the normal p + A suppression pattern 
is extrapolated to A + A That this is not theoretically fool proof was pointed 
out by Qiu et al|p8|. They showed that including radiative energy loss in cold 
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Fig. 14. Ratio of vr" yields in central to 
peripheral collisions at RHIC reported 
by PHENIX H scaled by the number 
of binary collisions Ncoii ~ 857, 19 resp. 
The suppression of high px pions at 
RHIC (in contrast to the Cronin en- 
hancement at SPS in Fig.|ll|) is due to 
jet quenching in the ultra-dense matter 
formed at RHICfRli^jTSHFlsiir 



pT[ai!V/cj 

Fig. 13. Semi-inclusive tt" pt distri- 
bution {l/Nint){dN^a /2TTpTdpTdy)) in 
the upper 60-80% peripheral events 
(solid circles) and the 10% most 
central events (solid squares) from 
PHENIX The lines are a param- 
eterization of pp charged hadron spec- 
tra, scaled by the mean number of col- 
lisions Ncou = 857, 19 resp. The bands 
indicate the possible range due to the 
systematic error on NcoU- 



nuclei could lead to non-linear enhancement of J/?/; suppression by decreas- 
ing their formation probability. Only a rather schematic model was presented, 
but it emphasizes the necessity of improving considerably the theory of the 
"normal" processes associated with heavy quark propagation through with 
nuclei. The situation is rather similar theoretically with regard to the Cronin 
effect. There also only rather schematic models are available to simulate the 
effect. 

The big difference between the two cases is that for J /if) the "normal" 
and "anomalous" components work in the same direction. The premium is 
thus high on developing an accurate theory "normal" nuclear suppression. 
In the jet quenching case, on the other hand, the "normal" Cronin effect 
works in the opposite direction to the "anomalous" new jet quenching mech- 
anism. Of course, there are possibly other "normal" effects, such as gluon 
(anti?)shadowing, that may work in either direction at high pT- To map out 
all the "normal" physics components will require detailed systematic mea- 
surements of p-f A at RHIC as done at the SPS. As a final remark, I want to 
emphasize the "normal" component of the dynamics is not dull run-of-the- 
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mill background, but fundamentally interesting many body QCD physics in 
its own right and deserves considerable more attention. 

3.4 Where Have all the Baryons Gone? 

One of the puzzling feature of Figs.(|l2|,^4|) is that pions appear to be more 
quenched than the sum of charged particles. Usually we assume that pions 
are the most abundant hadron species at high pT since both quark and gluon 
fragmentation functions prefer to make the lightest mesons [^J92| . Surpris- 
ingly, the preliminary PHENIX data[Q on identified high pT hadron spectra 
suggest from Fig.|l^ that baryons may be the most abundant species above 
Pt > 2 GeV/c. One possible source of such a non-pQCD like flavor distribu- 
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Fig. 15. Minimum bias transverse momentum distributions for positive (left) and 
negative (right) identified hadrons measured in PHENIX|Q. The error bars include 
statistical errors and systematic errors in the acceptance and decay corrections. 
Additional 20% systematic errors on the absolute normalization are not included. 



tion could be hydrodynamic transverse flow. For a longitudinal boost invari- 
ant (Bjorken) expansion with a transverse flow velocity field, v±_ = ta'Dh{rir), 
the general formula [ p3| for the differential particle number is 

dNs d J drf f mT(Mr\ f PTshr]r\ 



18 Miklos Gyulassy 



+^PtKo (^) (^) I , (6) 

where d = 2s + 1 is the degeneracy factor, rjj. = Artanh(v^(z = 0)) is the 
transverse fluid rapidity and (''/(C)j '''/(O) is a parameterization (counter- 
clockwise) of the freeze-out surface (isotherm of temperature Tf). 

Solutions for freeze-out surfaces with arbitrary transverse velocity fields 
f_L(0 '^^^ be obtained by solving relativistic hydrodynamics. For the simplest 
case with vj_ = tanhyy^ a constant and an isotherm freeze-out on a proper 
time hypersurface r/ Q, 

dNs d 2 ( '^T cosh r\r \ r f ^ii^^ "Hr 

;R TflTlTKl Iq 



dyd^PT 47r2 ^ V Tf J ' \ Tf 

PT»m ^^^^^ ^ ^ / _ PT ^ \ j--7\ 

V Tfexp{i]r)J 

which corresponds to a blue shifted effective temperature TfC^''. This is the 
uniform rapidity, transverse boosted Bjorken sausage parameterization of nu- 
clear collision distributions. 

Evidence for increased transverse flow phenomena at RHIC relative to 
SPS comes from low pT STAR data|9| shown in Figs.(|l|,0). The data can 
be fit up to < 1 GeV/c with a rather radial flow velocity v± ~ 0.6 c that 
is significantly larger than the radial flow ^ 0.4 c deduced from similar SPS 
spectra. 

Another important experimental tool to search for collective flow effects 
is to study anisotropic multiparticle emission patterns | p5| , p6p7|j9^ . A par- 
ticularly useful measure of collective behavior in ultra-relativistic energies 
has turned out to be the differential second Fourier component |Q of the 
azimuthal distribution: 

dNh(h) dNh(h)l , „ , 

where (j) is measured relative to a "reaction plane" event by event as deter- 
mined in Jq^ . 

Azimuthal or "elliptic" flow results from the initial spatial anisotropy 
of the dense matter in semi-peripheral A + A collisions. The hydrodynamic 
model predicts an elliptic flow pattern at RHIC |5^,|5^ that can be approxi- 
mately parameterized as 

w^(pT)«tanh(pT/(10±2GeV)) . (9) 

Up to about pt < 1 GeV, this agrees remarkably well with STAR data. 
At high pt this hydrodynamic flow component breaks down because of the 
emergence of the hard pQCD hadrons. 

The transverse boosted Bjorken sausage model eq.(|^) predicts that asymp- 
totically the baryon/meson ratios p/tt-\- = p/tt~ — s- 2 for any flow velocity 
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Fig. 16. Transverse momentum dis- 
tribution for n,K,p measured by 
STArQ at RHIC . The dashed curves 
are fits with eq.^ leading to a freeze-out 
temperature T/ = 100 MeV boosted 
with a radial flow velocity fleld v_i_ ~ 
0.6c. 




il 0.3 0.4 0.5 U 0.7 0.8 0.9 1 
p,[GeV/c] 

Fig. 17. Second azimuthal Fourier 
component of invariant pion and pro- 
ton distributions from STAR[^ com- 
pared to ideal hydrodynamic flow pat- 
tern from ISql. 



because nucleons have 2 spin states. However, numerically this ratios exceed 
unity only > 3 — 4 GeV. Thus transverse flow alone cannot account for 
the anomalous baryon dominance of high px spectra in Fig.|l^ as emphasized 
in @. 

Another observation |Q that possibly provides a hint that the answer to 
the puzzling result may lie in novel baryon dynamics at RHIC can be seen 



in Figs.( |l8| , pi)D . As was shown by Kharzeev|101|, the energy and rapidity de- 
pendence of the inclusive baryon production at mid-rapidity can be obtained 
using Mueller's generalized optical theorem in the double Regge limit. Here, 
the exchanges of a Pomeron and a Mq Baryon-anti-baryon "junction" pair 
lead to the following form for single mid-rapidity baryon production, 

EB-^=CBfB{fn't)(^) cosh(y/2) . (10) 



where Cb is a constant that reflects the couplings of the Reggeon and Pomeron 
to the proton, /^(toj) is an unknown function of and sq — 1 GeV is a 
Regge energy scale. The cosh(j//2) rapidity dependence and the 1/ energy 
dependence follow from the assumed intercept [ |102| , a^jj(O) ~ 1/2. In con- 
trast to simpler diquark breaking models as in the Dual Parton Model, the 
multiplicity of junction also enhanced events is enhanced by a factor of 5/4 
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Fig. 18. The valence proton rapid- 
ity density measured by STAR[^ at 
RHIC as a function of Npart- prelimi- 
nary BRAHMS data are also indicated. 




Fig. 19. Predicted valence proton ra- 
pidity density at RHIC from ]lOO 
showing factor of two enhancement ex- 
pected if baryon junction exchange is 
included in HIJING/B. The dashed 
curves are result of HIJING including 
only standard LUND diquark fragmen- 
tation. 



m.p + p, and the strangeness content is also enhanced by a large factor. The 
junction mechanism for baryon number (vs valence quark number) trans- 
port predicts for the unique possibility of producing S — fi~ baryons at 
midrapidity, as were observed at the SPS in WA97. In the Monte Carlo event 
generator HIJING/ BB ^5 10C|, baryon junctions are implemented in terms 
of Y shaped strings spanning valence quarks. 

The junction is a topological knot in the gluon field connecting the color 
flux from three quarks into a color singlet state[102|. The intriguing aspect of 
junctions is that the cons erve d baryon number resides in the gluon knot and 
not in the valence quarks|101|. In a nuclear collision some or all of the valence 
quarks may fragment into mesons. However, the gluonic junctions insure that 
baryon number is conserved. The understanding of the dynamics of junction 
exchange and pair production is still rather primitive, but the consistency of 
the baryon stopping power at SPS and now RHIC with HIJING/ BB predic- 
tions suggest that baryon dynamics at central rapidities may be especially 
interesting at RHIC. See ref. |102] for a discussion of possible novel junction 
network physics that may lead to femto-scale buckyball and even CP odd 
junction network production in A + A. 
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3.5 Quenching of Elliptic Flow 

As seen in Figs.(|l^,^) strong elliptic flow was discovery at RHIC consistent 
with hydrodynamic prediction at low transverse momentum px < 1 GeV. 
However, the preliminary data from STAR p04| shows that above pT > 2 
GeV the elliptic flow saturates and the azimuthal asymmetry deviates more 
and more from hydrodynamic behavior as seen in Fig.El^. This information 
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partons in a gluon plasma with rapidity 
density dNg/dy = 200,500,1000 from 
0- 




provides insight into how hydrodynamic behavior breaks down at high px due 
to the finite energy loss of partons in the plasma. As shown in detail in 
the saturation pattern at high px depends on the energy dependence of the 
gluon energy loss as well as on the geometry of the plasma density at finite 
impact parameters. It therefore provides tomographic information about the 
density profile and its evolution in ^ -|- ^. See sec 4.4 for more details. 



3.6 Where did the slowly burning plasma log vanish? 

The last major RHIC result that I highlight here is on pion interferometry. 
Relativistic combustion theory] 105, 10^, 93| predicts that if there were a suffi- 
ciently rapid cross over between the QGP and hadronic phases of ultra-dense 
matter, then a deflagration burn front may appear between two phases. The 
main characteristic of that burn front is its very small velocity in case the 
entropy density jump across it is sufficiently large and no high degree of non- 
equilibrium supercooling arises. Even with a smooth cross-over transition, 
such slowly burning plasma solutions were shown to exists as long as the 
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width of the transition region is {ATc/Tc < 0.08). The hfetime of a Bjorken 
plasma log is therefore significantly enhanced r ~ R/vd, where Vd ^ 1/25 is 
the small deflagration velocity in the static I+ID case. 

This characteristic time delay of the hadronization from a QGP state 
was suggested in 



107,1 



to be testable via pion interferometry. In ref. [|93[ 
the 3+lD hydrodynamic equations were solved to study this plasma "stall" 
phenomenon in detail. 

The two pion correlation function measures the coincidence probabil- 
ity P(pi,p2) of two (identical) bosons with momenta pi, p2 relative to the 
probability of detecting uncorrelated particles from different events, 



C2(Pl,P2) 



^(Pl,P2) 
^'(P1)^(P2) 



(11) 



If the average 4-momentum is denoted as K^^ — (p^ +p2)/2 and the relative 
4-momentum as = Pi ~ P2 j then under the assumption that the particle 
source is chaotic and sufficiently large. 



C2(Pl,P2) = 1 



■J^ dS-K exp[iIJ-q] f{^) 



dN 



^ d3p2 



where 



E 



dN 



d3p (27r)3 



dS -p 



(12) 



(13) 



is the single inclusive momentum distribution, f{x) = (e^ — 1)""'^, and u'^ the 
fluid 4- velocity. The integrals run over the assumed freeze-out hypersurface. 
In general, that hypersurface is represented by a 3-parametric (4-vector) 
function -£'^(C, Vi 4')i the normal vector on the hypersurface is determined 

by 



dE„ 



dS'^ dSf^ dS^ 
d( dr] dcf) 



dC drj ( 



(14) 



where tfiafi-y is the completely antisymmetric 4 tensor. For the common 
isotherm freeze-out temperature Tf hypersurface, the fluid velocity gerenally 
varies = u'^{IJ). 

For the Bjorken cylinder geometry, it is useful to restrict consideration to 
particles emitted at midrapidity, = — 0. Rotational symmetry around 
the z-axis in central collisions makes it possible to choose the average trans- 
verse momentum as K.± = {K, 0, 0), and consequently, C2{K, gout, gsidc) is a 
function of three independent variables only. The so called out and side pro- 
jections of the relative momenta are qout = (9out,0, 0), qsidc — (0, gsido,0). 
As shown in |107 10S| the width, l/Rgide, of the correlation function in gsidc is 
a measure of the transverse decoupling or freeze-out radius, while the width 
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1/Rout of the gout correlation function is also sensitive to the duration of 
hadronization, At 

Thus a QGP stall would manifest experimentally in Rout ^ Rside- In it 
was found that for possibly realistic parameters, Rout I Rside ~ 2 — 3, could 
be observed if a QGP stall occurred. 
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Fig. 22. Preliminary PHENIX pion interferometry data ]109| vs different projec- 
tions of the relative momenta. Similar preliminary STAR data|^ were also shown 
at QMOl. Unlike predictions^] there is no hint of the expected stall or time delay 
of the QGP transition. 



With this "warm-up" review of pion interferometry, we are now ready for 
the rude awakening from ideal gedanken considerations with the first splash 
of "cold" RHIC pion interferometry data shown in Fig. ( p2| ) . The preliminary 
PHENIX data show that Rout ~ Rside and even more disturbing all the 
deduced interferometry parameters are virtually identical to values seen at 
the AGS and SPS. To add insult to injury, it appears that Rout < Rside for 
Pt > 0.4 GeV. Preliminary STAR dataj^ show the same tendency. 

Of course scenarios may be invented to "explain" the data a postiori, but 
if these data are confirmed by further measurements, then they are indeed 
surprising and call into question our picture of the space-time evolution of 
A + A. That this problem is not restricted to idealized hydrodynamics can 
be seen from the results of Ref. |110|. It was shown that Rout > Rside is also 
predicted in a calculation where the entropy jump is small and pion decou- 
pling is dynamically handled via URQMD. Among the theoretical questions 
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that should now be further investigated is whether the pion interferomctry 
theory based on chaotic ensembles|lll] is in fact apphcable to A + A. An- 
other question that needs further study is whether the assumed ensemble of 
initial conditions was too restrictive and whether highly inhomogeneous and 
turbulent initial conditions apply pO]. 



4 Jet Quenching: Theory 

Having had a brief tour of some of the interesting new data harvested from 
RHIC during the first round of experiments, I turn next to the more specific 
theoretical problem of computing the energy loss per unit length of a fast 
parton penetrating a finite, expanding quark-gluon plasma. As I emphasized 
above, high pt many body pQCD physics is a new frontier at RHIC and 
higher energies. This requires the development of the non-abelian analogue 
of the radiative energy loss theory familiar from classical E&M. The inter- 
esting new twist is that we have no external beams of quarks or gluons and 
the medium is very thin due to the fact that nuclei are tiny. Also the for- 
mation time physics of Landau-Pomeronchuk-Migdal (LPM) results in major 
destructive interference effects that must be taken into account. Work on this 
problem over the past five years has advanced considerably but many open 
problems remain. 

I will only highlight only one of those direction, namely the opacity ex- 
pansion reaction operator method that we developed in refs.]?^. The reader 
is referred to BDMSQ, Z||7|l, and U(7|] for alternative methods and ap- 
proximations. 

In ref . Jz^ we proposed a simple model to study induced gluon radiation 
due to multiple elastic scattering of a high energy jet propagating in a locally 
color neutral amorphous plasma. 



4.1 GLV Formalism 

In ifzol we developed a systematic graphical method to compute medium in- 
duced gluon radiation amplitudes as shown in Fig.|2^. The exponential growth 
of the number of graphs with the number of interactions makes it very te- 
dious to go beyond order three. In GLV ||75| we overcame the combinatorial 
problem by developing a new algebraic operator technique to solve the inclu- 
sive radiated gluon distributions recursively. The first step is to compute the 
three direct (single Born) and four surviving virtual (contact double Born) 
diagrams shown in Fig|2J. 

For scattering off of n scattering centers located at depths Zi in a trans- 
verse homogeneous medium of large area ((/ii?)^ ^ 1), we can write the 
inclusive radiated gluon spectrum, P„ (k, c) , as a sum over products of par- 
tial sums of amplitudes and complementary complex conjugate amplitudes. 
Every term in the sum contributes to the same 0{g^^). The average value of 
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Fig. 23. Induced radiation amplitude [[7^] contributing to fifth order and higher 
order in the opacity expansion of QCD energy loss in the GW model 111. The 
crosses denote static color screened Yukawa interactions on a scale /x. The blob is 
the initial hard jet amplitude without final state interactions. 




Fig. 24. Three first order (singe Born) direct and four surviving (double Born) 
virtual or contact amplitudes |fr7| from which the Dn and Vn components of the 
reaction operator in eq.(p^) are derived in 
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Fig. 25. Example of graphs constructed via Di, Vi that contribute to the 4th order 
in opacity in elastic and inclusive inelastic final state interactions. The longitudinal 
depth of active scattering centers are denoted by Zi and inactive (created with li) by 
{zi). The form of Dt, Vi depend on the process type but the tensorial bookkeeping 
of partial sums of amplitudes is the same. 



n is referred to as the opacity of the medium. The partial sums of diagrams 
at order n in such and opacity expansion can be conveniently expressed in 
a tensor notation and constructed by repeated operations of, or Vi 

corresponding to no, direct, or virtual interactions at scattering center i 

n 

Ai - i„(a;,k, c) = Y\_ (^o,i„ +52,z„V:m) Go(x,k, c) . 

m— 1 

Here Gq is the initial hard q + g color matrix amplitude. In the inclusive 
probability each class contracts with a unique complementary class 

with the complementary class constructed as 

n 
m=l 



Lectures on the Theory of High Energy A+A at RHIC 



27 



Fig. ( p5[ ) shows an example of how this formahsm works at 4th order in opacity 
for elastic and inelastic inclusive distributions. 

Direct interactions enlarge rank n — 1 class elements as follows: 

-^n^zi k, c) — [cLji + Sn + Bji^jAi-^ ...i^ -^ (x, k, c) 

— Q^n-^ii •••in — i (*^; k, c) -t- e ^ ^ ^ -^ii — \ ^ i [^7 ^n] ) 

[ -2 j B„e*"°""[c,a„]Tei(Ai...«„_J 

where B„ = H — C„ = k/k^ ~ (k — qri)/(k — q™)^) is the so-called Bertsch- 
Gunion amplitude for producing a gluon with transverse momentum k in 
an isolated single collision with scattering center n. The momentum transfer 
to the jet is q„. The notation a;„ = (k — q„)^/2a;, for a gluon with energy 
w and a„ is the color matrix in the dn dimensional representation of the 
jet with color Casimir Cj^. = J2?n^i^im,2 counts the number of virtual 
interactions in Aij^...i^_^. 

Unitarity (virtual forward scattering) corrections to the direct processes 
involve the sum of four double born contact diagrams in Fig. ( p^ that enlarge 
rank n — 1 classes as follows: 

Vn = ~ 2^^^ + ~ ^ anBn = -a„Ai - 2^(^^ ~ ^^^^ 

This key operator relation between direct and virtual insertions that we dis- 
covered in 1^ makes it possible to solve the problem algebraically. 

The tensor classification of classes of diagrams makes it possible to con- 
struct the distribution of radiated gluons in the case of n interactions, P„, 
recursively from lower rank (opacity) classes via a "reaction" operator 

P„ = A''-'"-'RnA,...,„_, , Rn = DiDn + K + (16) 
Using the key identity (^5|), the reaction matrix simplifies to 

Rn = (A. - an)HDn - a„) - Ca^ [Sn + B„)U^„ + S„) - Ca 

The next major simplification occurs because both S and B involve the 
same gluon color rotation through This reduces the color algebra to 

multiplicative Casimir factors 

= Ca (P„-i(k - q„) - P„-i(k)) = Ca [e'^" '' - l) P„-i(k) 
Jn - i„-i Rt n 4. — n 
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In obeys a recursion relation from which the inclusive radiation probability 
is found to obey the soluble recursion relation 

P„(k) =C^(P„_i(k-q„)-P„-i(k))-2CAB„. (Re e-'^-'-e'^" hn-i) 

+<5n,lC'ACfl|Bip 

where b = «Vk is the transverse momentum shift operator. The initial condi- 
tion for this recursion relation is the initial hard vertex radiation amplitude 
without final state interactions that is given by Pq — Cr = C/j/k^. 

The complete solution to the problem can therefore be expressed in 
closed form as 

n I n I 

P„(k) = -2CflQRe^ I Yl (e'"''^ - 1) \ • e'^^^-^e-^^^'^x 
i=i J 

This expression can be averaged over any spatial distribution of interaction 
centers, Zi as well as any Zi dependent momentum transfers g„. This form 
is thus ideally suitable for Monte Carlo implementation for arbitrary q_i,Zi 
medium ensemble averages. 



4.2 Nonabelian Energy Loss at Finite Opacity 

The first application of our general solution to the energy loss problem 
was to calculate numerically the total radiated energy loss as a function of jet 
energy E, plasma depth, L, and infrared screening scale, /i. In the absence of a 
medium, the gluon radiation associated with a spin ^ parton jet is distributed 
as 

,^^^(l-. + ^]± (17) 

where x — k'^/E^ w ^^/E, and Cu is the Casimir of the (spin 1/2) jet in 
the dji dimensional color representation. The differential energy distribution 
outside a cone defined by k^^ > is given by 

^ = ^^(,^-"+yj (18) 

where the upper kinematic limit is k^max = min [4£^^a;^, 4£'^x(l — x)] .. 
The energy loss outside the cone in the vacuum is then given by 

AE("^^^^E log- (19) 
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While this overestimates the radiative energy loss in the vacuum (self-quenching), 
it is important to note that /E ^ 50% is typically much larger than 

the medium induced energy loss. 

Averaging over the momentum transfer qi^ via the color Yukawa po- 
tential leads to a very simple first order opacity result for the a; ^ 1 gluon 
double differential distribution 

^-1 — TTT^^-i — ITT ^ "qi- 



dxdk^ dxdk^ Xg Jo n{qlj_+fi'^y 
2k^.q^(k-qi)i£^ 

16x^E^ + (k-(ii)]_L'^ ■ ^ ' 

where the opacity factor L/\g ~ Na'fi^ /A± arises from the sum over the N 
distinct targets. Note that the radiated gluon mean free path Xg = {Ca/Cr)X 
appears rather than the jet mean free path. The upper kinematic bound on 
the momentum transfer q^^^ = s/4 ~ SEfi, (l/^^^y = ^ +'7max))- 
For SPS and RHIC energies, this finite limit cannot be ignored as we show 
below. 

The second order contribution in opacity involving the sum of 7^ direct 



and 2 X 86 virtual and results in |75 



P(2^ cx CnCldR [ 2 Ci • Bi (1 - cos(tJiZ\zi)) 
+ 2 C2 • B2 {cos{uj2Az2) ~ cos{LU2iAzi + AZ2)) 

~ 2 C(i2) • B2 (cOs(Ci;2^^2) - C0S(W(12)Z\Z1 + LU2AZ2)) 

-2C(i2) •B2(12) (l-cos(w(i2)Z\zi))] , (21) 

where with C(„i„) and oj^rnn) obtained from H and cjq through the substitu- 
tion k_L ^ - q^,„ - q_L„ and B„(„;) = C,n ~ C(„;)[Q. 

Numerical results comparing the first three orders in opacity corrections 



to the hard distribution Eq. (^7|) were presented in |75 . To illustrate the 
result consider a quark jet in a medium with Ag = 1 fm, a screening scale 
/X = 0.5 GeV and = 0.3 . The total radiative energy loss could be expressed 
as 

/i^^«-§^^log- , (22) 

N{E) Xg 

with N{oo) = Alog{E/ fi)/v if the kinematic bounds were ignored as in the 
approximations of ref . JttI . We found that finite kinematic constraints and 
the form of the first order result cause N{E) to deviate considerably from 
the asymptotic value for all energies accessible in the RHIC range. Together 
with the logarithmic dependence on energy, these kinematic effects suppress 
greatly the energy loss at lower (SPS) energies as seen in Fig.|2^. This is in 
sharp contrast to the approximately energy independent result in BDMS-ZW 
where the finite kinematic bounds were neglected because only the asymptotic 
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Fig. 26. The GLV radiated energy lossQ of a quark jet witli energy Ejet = 
5, 50, 500 GeV (at SPS, RHIC, LHC) is plotted as a function of tiie opacity L/Xg. 
(Ag = 1 fm, fj, = 0.5 GeV). Solid curves show first order, while dashed curves show 
results up to second order in opacity. The asymptotic energy loss (solid triangles) 
of BDMS|^^ is shown for comparison. The energy dependence of GLV suppressing 
radiative energy loss of low energy jets explains why no jet quenching was observed 
at the SPS (see Fig.[o|) . 



limits were considered. Another remarkable result demonstrated numerically 
is that the second and third order contributions to the integrated energy loss 
remains surprisingly small in the physical range of nuclear opacities L/Xg ~ 5. 
The rapid convergence of the opacity expansion even for realistic opacities 
results from the fact that the effective expansion parameter is actually the 
product of the opacity and the gluon formation probability Ljj.^ /2xE.The 
leading quadratic dependence of the energy loss on nuclear thickness discov- 
ered in BDMS l??! therefore already emerges from the dominant first order 
term in the opacity expansion. 

At SPS energies kinematic effects suppress greatly the energy loss relative 
to BDMS. Our estimates provide a natural explanation for the absence of jet 
quenching in Pb + Pb at 160 AGeV observed by WA98. At RHIC energies, 
on the other hand, a significant nonlinear (in A) pattern of suppression of 
high p± hadrons relative to scaled pp data is predicted. 

4.3 The Opacity of the QGP at RHIC 

As a second a application of the GLV energy loss, in ref. we computed the 
quenched pQCD distribution of high px hadrons as a function of the effective 
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static plasma opacity, L/Xg. In Figs.( |27| , p8| ), the jet energy dependence of the 
GLV energy loss for gluons is shown. The most important feature to note is 
that AEqlv/E is approximately constant in the energy range accessible at 
RHIC. 
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The relative energy loss 
(AE/E) is approximately constant at 
medium energy, 2 < i? < 10 GeV. 



In order to compute the pion spectrum, note that jet quenching reduces 
the energy of the jet before fragmentation. We concentrate on mid-rapidity 
iUcm = 0), where the jet transverse momentum before fragmentation is shifted 
by the energy loss as in pT^ , p*iL/X) = Pc - AE{E,L). This shifts the 
Zc parameter in the fragmentation function of the integrand ( p3| ) to z* ~ 
2c/(l - AE/p,). 

The invariant cross section of hadron production in central A-l-A collision 
is then given bypl| 



E^ 



do^ 
d^p 



j dV iA(b)tB(b - r) ^ y"da;ada;bdzcd2fc±,ad2fci,b 

abed 



/a/A(a;a, fc±,a(b), Q^)h/A{xh, fcj.,b(b - r), Q2) 



da 
df 



s(5(s + t + u) , 



(23) 



where upper limit of the impact parameter integral is fcmax = 4.7 fm for 
10 % central Au-|-Au collisions. Here tj>,{h) is the usual (Glauber) thickness 
function. The factor z*/zc appears because of the in-medium modification 
of the fragmentation function [112|. Thus, the invariant cross section ( p3| ) 
depends on the average opacity or collision number, n = L/X^. The calculated 
spectra for pions are displayed for n — 0,1, 2, 3, 4 in Fig. Fig. ^ shows 
their ratios to the non-quenched spectra at fi = 0. We note that in contrast to 
previous energy independent estimates for the energy loss, the GLV energy- 
dependent energy loss leads to constant suppression of the high px domain in 
agreement with the preliminary data. The peripheral collisions are consistent 
with a rather small opacity in contrast to central collisions, as expected. 
The ratio of central to peripheral PHENIX [p5[ data from QMOl shown in 
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Fig. 29. Pion production in Au+Au 
collision including jet quenching with 
opacity L/A = 1,2,3,4. Preliminary 
QMOl PHENIX data shown (see up- 
dated data from ||5| in Fig.|l|) 
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Fig. 30. The ratio of the central to 
the peripheral pion yields (normalized 
by the number of binary collisions, 
857 and 5.5). (Note updated data are 
shown in Fig.[l4|) 



Fig. 30 clearly reveals that jet quenching at RHIC overcomes the Cronin 



enhancement at zero (final state) opacity. This is in stark contrast to data at 
SPS energies, where WA98 found no evidence for quenching in Pb + Pb at 
160 AGeV but a factor of two Cronin enhancement as discussed before. 

Figs. p9| , ^ indicate that an effective static plasma opacity i/A = 3 — 4 
is sufficient to reproduce the preliminary jet quenching pattern observed at 
RHIC. In ref. it was shown that a rather small constant dE/dx ~ 0.25 
GeV/fm was also found to be consistent with the data. However, it is impor- 
tant to emphasize that these effective static plasma opacities and parameters 
hide the underlying rapid dilution of the plasma due to expansion. The GLV 
formalism including the kinematic constraints at first order has been further 
generalized to include effects of expansion in ||2l[] . It was found in that 
the inclusion of longitudinal expansion modifies the static plasma results in 
such a way that the moderate static plasma opacity actually implies that the 
produced mini-jet plasma rapidity density may have reached dNg/dy ^ 500. 



4.4 Jet Tomography from Quenched Elliptic Flow 

So far we have not included the dilution effect of expansion on the energy 
loss. The generalization of GLV to the case of expanding plasmas is[^ 

OO |k|niax |q|max 

diGLV 9CrE , . f a, f d^^al 



J dz p{z) J d^kas J 



dx TT^ y J J (q^ -t- ^Ji.{zY) 
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kq 



k2(k-q; 



(24) 



where p{z) is the plasma density at time z along the jet path at position 
z from the production point and where the screening scale ^i{z) may also 
depend on time. 

Consider a density evolution of the form as in ||77| , 

p[z)^p,{^^ye{L-z), (25) 

where a = corresponds to a static uniform medium of thickness i, while a = 
1 to Bjorken 1+lD longitudinal expansion transverse to the jet propagation 
axis. 

Analytic expressions can only be obtained again for asymptotic jet en- 
ergies when the kinematic boundaries can be ignored. In that case, all but 
the path integral can be done giving 

« — — / dz —— {z - zo) log , 26 

2 A(z) Lp^{z) 

which is a linear weighed line integral over the local transport coefficient 
(/x^(z)/A(z)) Ri ^Tralp(z) however enhanced by a \og2E / L p? {z) factor that 
results from the structure of the GLV integral missing in the BDMS asymp- 
totic limit. For an expanding plasma as in ( p5| ) 

Cna, pHL)L- L^"" , 
AE4L,zo)^— ^^v. (27) 

Here v = log2E/ Lp^{L) and we used that p?{L)L°' / \{L) is a constant in- 
dependent of L for this type of expansion. We also took the zq ^ limit. 
We therefore recover the asymptotic BDMSZ energy loss for both static and 
expanding media modulated by a log£^/a;(L) factor that is important at 
RHIC energies. Using the Bjorken relation between the gluon density and 
the rapidity density then gives 

In practice, it is straight forward to integrate GLV numerically including the 
finite kinematic constraints. 

For non-central collisions the GLV line integral depends of course on the 
azimuthal direction <f> of the jet. The variation of the azimuthal energy loss 
with respect to at a given impact parameter b can be expressed in terms of 

i.(b,0)^^f^'f) 

with results shown in Fig.O The effect of this azimuthal variation of the 
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Fig. 31. The modulation function _R(b, 0) is plotted vs. (p for impact parameters 
b = 2, 6, 10 fm. Diffuse Wood-Saxon and uniform sharp cylinder geometries are 
compared. The most drastic difference between these geometries occurs at high 
impact parameters. 



energy loss is to induce an apparent elliptic flow at high px not related to hy- 
drodynamic phenomena of low p^- In | pT| we proposed a simple interpolation 
between the hydrodynamic and jet quenched px eikonal regimes 



V2{pt) 



V2s{PT)dNs + V2h{PT)dNh 

dNs + dNh 



(29) 



This interpolates between the hydrodynamic and the pQCD regimes because 
at highpT, dNh 3> dNs- For our numerical estimates the lowpx interpolation 
is achieved by turning off the pQCD curves with a switch function fh{PT) = 
0.5[1 + tanh(3(pT - 1-5 GeV))]. 

We see in Fig.^ that the magnitude and shape of the high pT ellip- 
tic flow provides a complementary probe of the initial gluon density and is 
also sensitive to the geometrical distribution of the plasma. The saturated 
V2 increases systematically with increasing initial plasma density and thus 
provides an important complementary constraint on the maximum initial 
parton density produced in central b = collisions. The consistency of the 
quenched elliptic flow in non-central with the central quench pattern will be 
very important to test when the final data become available. 
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Fig. 32. The interpolation of V2{pt) between the soft hydrodynamic and hard 
pQCD regimes is shown for different gluon rapidity densities in central 6 = 
collisions. The gluon rapidity density at fe 7^ is assumed to scale with the binary 
collision number. Solid (dashed) curves correspond to sharp cylindrical (diffuse 
Wood-Saxon) geometry 



5 Summary 

If confirmed by further measurements and theoretical refinements, jet quench- 
ing may have already provided the first evidence that initial parton densities 
on the order of 100 times nuclear matter density may have been produced 
at RHIC. The full analysis of flavor composition, shape, and azimuthal mo- 
ments of the high pT spectra appears to be a promising diagnostic probe of 
the evolution of the gluon plasma produced at RHIC. However, it is too early 
to tell what the preliminary say about the properties of that extremely dense 
form of matter. There are too many pieces of the puzzle that simply do not 
fit well into any scenario. The beam energy and centrality independence of 
the transverse energy per charged particle is one of them. The anomalous 
baryon number transport to high transverse momenta and central rapidities 
is another. Finally, the puzzling beam energy independence of the prelim- 
inary pion interferometry results is a mystery. As the tera-bytes of RHIC 
data continue to stream in during the next few years, they will certainly pose 
many interesting new QCD many body problems. The new chapter on the 



36 Miklos Gyulassy 



physics of ultra-dense matter and the dynamics of uhra-relativistic nuclei is 
now unfolding at RHIC. 
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